SOME EXACT SOLUTIONS OF THE DIRAC EQUATION 
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Exact analytic solutions are found to the Dirac equation for a combination of 
Lorentz scalar and vector Coulombic potentials with additional non-Coulombic 
parts. An appropriate linear combination of Lorentz scalar and vector non- 
Coulombic potentials, with the scalar part dominating, can be chosen to give exact 
analytic Dirac wave functions. 



In a previous letter(!l simple exact solutions were found for the Dirac equa- 
tion for the combination of a Lorentz vector Coulomb potential with a linear 
confining potential that was a particular combination of Lorentz scalar and 
vector parts. In this work, we extend the method of Ref. 1 to the more general 
case of an arbitrary combination of Lorentz scalar and vector Coulombic poten- 
tials with a particular combination of Lorentz scalar and vector non-Coufombic 
potentials. A more complete version of this work will appear elsewhereo 

The Dirac equation we solve is 
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where a and (3 are the usual Dirac matrices. The four component wave function 
■0 can be written in terms of two component spinors u (upper component) and 
V (lower component) satisfying the equations 
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The key step in generating relatively simple exact solutions of the Dirac equa- 
tion is to choose a particularly simple form for the function v{r). For s-states,^ 
we choose v(r) = i"f{cr ■ f)u(r), where 7 is a constant factor to be determined 



"States with orbital angular momentum are treated in Ref. 2. 
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by the solution to the Dhac equation. This is the form of v{r) that was found 
in Rcf. 1 for a Coulomb plus linear confining potential. This ansatz for v{r) has 
also been used as the basis for generating approximate saddle point solutions 
for the DiracI and Breiu equations. Using this form of v{r), equations (2) and 
(3) reduce to two first order ordinary diffential equations for u{r) 
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Equations (4) and (5) are two independent equations for the same quantity, so 
that each term in one equation can be equated with the corresponding term 
in the other equation having the same radial dependence. This leads to 
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The relations in Eq. (6) can be rearranged, after some algebra, to give 
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The constant h can have either sign. Although b must be positive in the pure 
Coulombic case, we will see that a negative b is possible if the Lorentz scalar 
potential S{r) is more singular at the origin than 1/r. The bound state energy 
can be written as 
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The wave function u{r) can be found by solving Eq. (5) to give 
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where the constant a is given by a = j{m + E) = zt\/m'^ — E"^ . The constants 
7 and a can have either sign if S{r) approaches a constant or diverges as r 
becomes infinite. The integral in Eq. (^) can diverge for any finite r, as long as 
the product of a times the integral diverges in the positive sense. The integral 
can also diverge at the origin or as r —s- 00 as long as the quantity in square 
brackets in Eq. (0) remains negative. 
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The last equality in Eq. (g) shows that in order for this class of exact 
solutions to apply, the Lorentz vector and scalar non-Coulombic potentials 
must have the same radial dependence and opposite sign, with the vector 
potential being smaller than the scalar potential. As long as this constraint 
is satisfied, the results in Eqs. (0)-(||) represent a complete exact solution for 
the ground state wave function and energy of the Dirac Hamiltonian given in 
equation (1). 

Because of the constraints imposed on the potentials by Eq. (^), the energy 
can be written purely in terms of one set of potentials or the other. The 
ground state energy can be specified by the Coulombic constants A and r/, 
with the other equality in Eq. (^) serving as a constraint on the non-Coulombic 
potentials V{r) and S{r). Or the ground state energy could be specified by 
the ratio V{r)/S{r) of the non-Coulombic potentials with the other equahty 
serving as a constraint on A and 77. Although the possibility of this class of 
exact solutions is limited by the constraints on the potentials, this still permits 
a wide range of non-Coulombic potentials. 

We now consider conditions imposed on the potentials and the wave func- 
tion parameters by the physical requirements that the potentials be real and 
the wave function normalizable. We see from Eq. (|^) that 7 must be real, and 
then from Eq. (0) that b must be real. This requires the Coulombic potentials 
to satisfy the condition 1 — + ij'^ > 0. The reality of 7 restricts possible 
bound state energies to the range —m < E < m. Note that negative energies 
can occur, but E + m cannot be negative. Also, E cannot equal ±m, because 
this would make 7 = 0, leading to a constant, unormalizable wave function. 
This condition on E, along with Eq. (^), means that V{r) must always be less 
in magnitude than S{r). 

We discuss the remaining conditions on the parameters in terms of three 
sub-classes of solution: 

1. The "normal" class of solutions has b, 7, and a all positive. In this case, 
we see from Eq. (|^) that the Coulombic potentials must satisfy the further 
condition A -I- 77 > 0. 

2. This sub-case has b negative, with 7 and a still positive. The constant b 
can be negative if the product aS{r) is positively divergent at the origin faster 
than l/r. Then each of Eqs. (7)-(10) hold just as for positive b, and the wave 
function is still normalizable. The states with positive and negative b are not 
two ground states of the same Hamiltonian because the potentials cannot be 
the same for each state. That is, either the Coulombic potentials or the non- 
Coulombic potentials must change to be consistent with a negative b. Sub-case 
1 with positive b transforms smoothly into the pure Coulombic solution as the 
non-Coulombic potential tends to zero everywhere. But this is not true for sub- 
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case 2 with b negative. This sub-case requires the non-Coulombic potential to 
be dominant at the origin, and so has no corresponding pure Coulombic Umit. 

3. This sub-case has a negative 7 and a negative a, while b can have either 
sign, as discussed in sub-cases 1 and 2 above. A negative a is possible if the 
non-Coulombic potential diverges or approaches a constant as r 00, so that 
the integral in Eq. (|^) diverges faster than r at large r. Since a is negative, 
the potential S{r) must be negative at large r. Then all of Eqs. (7)-(10) hold 
as for positive a, and the wave function is still normalizable. This case is 
highly unusual, because it allows the possibility of a potential that is negative 
everywhere and diverges negatively at both the origin and infinite r. We know 
of no other example in quantum mechanics where such a negative potential 
can lead to a normalizable ground state. The reason this is possible here can 
be seen from Eqs. (4) and (5). There it is seen that S{r) enters the differential 
equations for u{r) only in the combinations "/S(r) or S{r)/"f. Since these 
effective potentials are positive, the resulting wave function is normalizable. 
As with sub-case 2, the case with 7 and a negative does not approach a pure 
Coulombic case if the non-Coulombic potential tends to zero. 

We now look at some special cases. If the non-Coulombic potentials are 
absent, then the solutions are for a general linear combination of Lorentz vector 
and Lorentz scalar Coulombic potentials. If either constant, A or 77, is zero, 
we recover the usual solutions of the Dirac equation for a pure scalar or vector 
Coulombic potential. The Coulombic potentials cannot both be absent (while 
keeping a non-Coulombic part) because then 7 would be zero leading to a 
constant, unnormalizable wave function. 

The method we have described does not work for radially excited states, 
because the simple ansatz for v{r) does not lead to consistent equations for 
du/dr in that case. The method does work for the lowest orbitally excited 
state for which I = j — ^. That case is discussed in Ref. 2. 
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